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ABSTRACT 

In this work we investigate the detectability of the gravitational stochastic background 
produced by cosmological sources in scenarios of structure formation. The calculation 
is performed in the framework of hierarchical structure formation using a Press- 

OSchechter-likc formalism. The model considers the coalescences of three kind of 
binary systems, namely double neutron stars (NS-NS), the neutron star-black hole 
^.J binaries (NS-BH), and the black hole-black hole systems (BH-BH). We also included 

in the model the core-collapse supernovae leaving black holes as compact remnants. 
In particular, we use two different dark-energy scenarios, specifically cosmological 
I constant (A) and Chaplygin gas, in order to verify their influence on the cosmic star 

formation rate, the coalescence rates, and on the gravitational wave backgrounds. We 
calculate the gravitational wave signals separately for each kind of source as well as we 
determine their collective contribution for the stochastic background of gravitational 
waves. Concerning to the compact binary systems, we verify that these sources 
produce stochastic backgrounds with signal-to-noise ratios (S/N) ~ 1.5 (~ 0.90) 
for NS-NS, - 0.50 (~ 0.30) for NS-BH, - 0.20 (~ 0.10) for BH-BH, and - 0.14 
(~ 0.07) for core-collapse supernovae for a pair of advanced LIGO detectors in the 
cosmological constant (Chaplygin gas) cosmology. Particularly, the sensitivity of the 
future third generation of detectors as, for example, the Einstein Telescope (ET), in the 
triangular configuration, could increase the present signal-to-noise ratios by a high- 
factor (~ 300 — 1000) when compared to the (S/N) calculated for advanced LIGO 
detectors. As an example, the collective contribution of these sources can produce 
(S/N) ~ 3.3 (~ 1.8) for the A (Chaplygin gas) cosmology for a pair of advanced LIGO 
interferometers and within the frequency range ~ 10 Hz — 1.5 kHz. Considering ET we 
have (S/N) - 2200 (~ 1300) for the A (Chaplygin gas) cosmology. Thus, the third 
generation of gravitational wave detectors could be used to reconstruct the history of 
star formation in the Universe as well as for contributing with the characterization of 
the dark energy, for example, identifying if there is evidence for the evolution of the 
dark energy equation-of-state parameter w(a). 

Key words: binaries: close - gravitational waves - black hole - large-scale structure 
of Universe - stars: neutron - cosmology: theory - dark energy. 
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1 INTRODUCTION 

The direct detection of gravitational waves (GWs) is a major 
challenge for the physics/astrophysics and considerable 
experimental effort is being devoted by several groups 
around the world. In particular, the window in the frequency 
range ~ 10 — 10 kHz is open due to pioneerism of the 



following interferometers: LIGO detectors (e.g., Abbott et 
|al.|2009| >, VIRGO detector (e.g. |Acernese et aL|2008| >, GEO 
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600 detector (e.g., |Grote et al.||2008[ ), and TAMA 300 
detector (e.g., [Takahash i~et al.|2004f l 

In the future, another window in the low frequency 
range ~ 10~ 4 — 10 Hz will be open by space antennas 
such as e-LISA/NGO ( jAmaro-Seoane et al. ||2012|), BBO 
( Cutler fc Harms|2006| ), DECIGO ( |Ando et al 



[20121 , 

poiof . 



These 

interferometers together ground-based detectors which are 
presently in project 
and 



the third 



as KAGRA (Somiya et al 
generation of 



2012), 



as SCHENBERG (Aguiar et al 



resonant mass detectors 
20081 and MiniGRAIL 
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( Gottardi et al.|2007 1 will transform the research in general 
relativity into an observational/theoretical study. 

On the other hand, as GWs are produced by a 
large variety of astrophysical sources and cosmological 
phenomena, it is quite probable that the Universe is 
pervaded by a background of such waves. Collapse of 
Population II and III stars, phase transitions in the early 
Universe, cosmic strings, and a variety of binary stars are 
some examples of sources that could produce such a putative 
background of GWs (see, e.g., |Maggiore|[2000 Rcgimbau 



& de Freitas Pachcco 2001; de Araujo, Miranda & Aguiarl 


2000[ |2002[ |2004 


dc Araujo <fc Miranda 2005; Miranda, de 


Araujo & Aguiar 


2004[ Sandick et al. 


2006 


Suwa et al.|2007| 


Giovannini 2009 


Pereira & Miranda 


2010 


among others). 



Observe that the indirect evidence for the existence 
of gravitational waves came first from observations of the 



orbital decay of the Hulse- Taylor binary pulsar (Hulse & 
Taylor 1974 1975a|b I. In this century, direct detection 
though and analysis of gravitational-wave sources are 
expected to provide a unique insight to one of the least 



understood of the fundamental forces (Belczyhski, Kalogera 
|fc Bul ik 2002). Specifically, gravitational waves could also 
be used as a tool for studying the viability of different 
alternative theories of gravity. This could be done by 
comparing the detected polarization modes with those 



predicted by general relativity (see, e.g., de Paula, Miranda 
fc Marinho||2004| |Alves, Miranda fc de Araujo|2009[ ). 

As mentioned above, a number of interferometers 
designed for gravitational wave detection are currently in 
operation, being developed, or planned. In particular, the 
high frequency part of the gravitational wave spectrum 
(10 Hz < / < 10 kHz) is open today through the pioneering 
efforts of the first-generation ground-based interferometers 
such as LIGO. While detections from this first generation 
of detectors are likely to be rare, the third generation of 
gravitational wave detectors as, for example, the Einstein 
Telescope (ET) may detect, among others, the stochastic 
signal generated by a population of pre-galactic stars. Thus, 
gravitational wave observations could add a new dimension 
to our ability to observe and understand the Universe. 

On the other hand, the state of the art in cosmology 
has led to the following distribution of the energy densities 
of the Universe: 4% for baryonic matter, 23% for non- 
baryonic dark matter and 73% for the so-called dark energy 
( Jarosik et al.pbll I. Concerning to the dark energy, some 
equations of state have been proposed in order to explain 
such a dark component. The most common example is the 
cosmological constant (ACDM model), which implies on a 
constant vacuum energy density along the history of the 
Universe. Another possibility is a dynamical vacuum or 
quintessence. In general, the quintessence models involve 



one ( Albrecht & Skordis 2000 \ or two ( Bento, Bertolami 
|fc Santos||2002| > coupled scalar fields. The Chaplygin gas 
is another example of dark energy fluid. One of the most 
appealing aspects of the original Chaplygin gas model is that 
it is equivalent to the Dirac-Born-Infeld (DBI) description 



of a Nambu-Goto membrane 1 


Bordemann & Hoppe 1994 


Gorini et al.|2005 Ogawa|2010 


1. 



In this way, the main goal of the present work is 
to explore the possibility of using stochastic backgrounds 
of gravitational waves to provide more information about 
the character and interrelationship of the dark-energy 



equation of state with the star formation at high redshifts. 
That is, in first place we analyze the influence of 
two different dark-energy components of the Universe, 
namely cosmological constant and Chaplygin gas, on the 
stochastic backgrounds of gravitational waves produced by 
four different cosmological sources which are the merging 
together of two neutron stars (NS-NS), the coalescence of 
neutron star-black hole (NS-BH) systems, the merger of 
two black holes (BH-BH), and the core-collapse supernovae 
leaving black holes as compact remnants. In second place, 
we show that different dark-energy fluids produce distinct 
signatures for the cosmic star formation rate (CSFR) 
specially at high redshifts (z > 3). This interesting feature 
could be used as an alternative way to study the star 
formation and the possible temporal dependence of the dark- 
energy equation of state up to redshift 20, having as the 
common tool a stochastic background of gravitational waves 
with high signal-to-noise ratio. In this way, not only binary 
systems at lower redshifts (z < 2 — 3) working as standard 
sirens but also stochastic backgrounds of gravitational 
waves could contribute for a better comprehension of the 
physical nature of the dark energy and their connection, 
and influence, with the star formation at high redshifts. 

The preference for concentrating attention on the 
Chaplygin gas also comes from recent work of |Pace,| 



Waizmann & Bartelmann ( 2010 \ who analyzed the spherical 



collapse model in dark-energy cosmologies. As can be seen 
from that work, the Chaplygin gas exhibits an equation- 
of-state dependent on time as the quintessence models also 
exhibit. It is not the purpose of the present study make an 
individual assessment of each particular type of dark-energy 
candidate. Our goal is to verify if stochastic backgrounds of 
gravitational waves can give us some indication about the 
evolution of the dark-energy equation of state with time. 
Thus, the comparison between Chaplygin gas and A-CDM 
is sufficient for the purposes of this study. 

Here, we start with the CSFR recently derived by 



Pereira & Miranda (20101. Specifically, these authors use 
a hierarchical structure formation model and they obtain 
the CSFR in a self-consistent way. That means, the authors 
solve the equation governing the total gas density taking 
into account the baryon accretion rate, treated as an infall 
term, and the lifetime of the stars formed in the dark matter 
haloes. Here, we adapted the formalism derived by [Pereira] 
fc Miranda] ( |2010[ ) in order to obtain the CSFR and the 
coalescence rates in consistency with the assumed dark- 
energy model. 

The paper is organized as follows. In Section 2, we 
review the basics of the hierarchical model and how to 
obtain the CSFR up to redshift z ~ 20 as a function of 
the specific dark-energy cosmology. In Section 3, we discuss 
how to obtain the coalescence rates for NS-NS, NS-BH, and 
BH-BH systems from the CSFR. In Section 4, we present 
the formalism used to characterize the gravitational wave 
backgrounds for compact binary systems and core-collapse 
supernovae to form black holes. We also present the signal- 
to-noise ratios (S/N) for both a pair of advanced LIGO 
detectors and the ET in triangular configuration. Section 
5 presents the collective contribution of these sources for 
the gravitational wave background. In Section 6 we discuss 
the influence of the uncertainties of the parameters on the 
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derived gravitational wave background and on the CSFR. 
Section 7 presents the final considerations of this work. 



2 THE COSMIC STAR FORMATION RATE 
AND THE DARK ENERGY COSMOLOGIES 

2.1 An overview 



The essence of the halo model was discussed by |Neyman^ 
Scott & Shane ( 1952 I who postulated that all galaxies form 



in clusters, the distribution of galaxies within clusters can 
be described by a probabilistic relation, and cluster centers 
are themselves correlated. Substituing the word 'clusters' by 
'haloes' in the paper of Neyman, Scott & Shane ( 1952 1 we 



arrive at a reasonable qualitative description of the modern 
halo model. Today, it is widely believed that haloes, or 
overdense dark matter clumps, form as result of the growth 
and non-linear evolution of density perturbations produced 
in the early Universe (Peacock 19991. This is the heart of 



the hierarchical formation scenario. 

In general, the halo mass function is represented as the 
differential number density of haloes with mass between 
M and M + dM. Press and Schechter (hereafter PS) 
heuristically derived a mass function for bound virialized 



objects in 1974 (Press & Schechter 19741. The basic idea 



of the PS approach is define haloes as concentrations of 
mass that have already left the linear regime by crossing 
the threshold S c for non-linear collapse. Given a power 
spectrum and a window function, it should then be relatively 
straightforward to calculate the halo mass function as a 
function of the mass and redshift. 

However, it is worth stressing that the exact definition 
of the mass function, e.g., integrated versus differential 
form or count versus number density, varies widely in the 
literature. To characterize different fits, it can be introduced 
the scale differential mass function 



(Jenkis et al. 

2001) defined as a fraction of the total mass per In a~ that 



belongs to haloes. That is, 



diner" 



M 



dn(M, 



P b(z) dMff-^M,*)]' 



(1) 



where n(M, z) is the number density of haloes with mass 
M, and Pb(z) is the background density at redshift z. As 



pointed out by Jenkis et al. (20011, this definition of the 



mass function has the advantage that it does not explicitly 
depend on redshift, power spectrum, or cosmology; all of 
these are contained in a{M,z) (see also Lukic et al.||2007 |. 
See that 



AC*) = 



J^ : f(a)MdM 
f °° f(a)MdM 



(3) 



where we have used M min = 10 6 M B and Af max = 10 M 



(see Pereira fc Miranda|2010 for details) 



Therefore, the baryon accretion rate cib(i) which 
accounts for the increase in the fraction of baryons in 
structures is given by 



db(£) = fibPc 





df b {z) 




dz 



(4) 



where p c — 3Hq /8ttG is the critical density of the Universe. 

The age of the Universe that appears in |4| is related 
to the redshift by: 



dt 
dz 



(l + z)E(z)- 



(5) 



In Eq. |5 |) E(z) represents the expansion function which 
is (see, e.g., Pace, Waizmann fc Bartelmann||2010 \ 



n m (l + z) 3 + Ojexp ( -:! / 1 + w ( ffl, ) rfQ 



E(z) 



(6) 

where the relative density of the i-component is given by 
Sli = pi/pc, and 'i' appying for baryons (b), dark energy 
(d), and total matter (m). As usual, the scale factor is 
a = 1/(1 + z), and w(a) is the dark energy equation-of-state 
parameter. 

Note that for w(a) = — 1 we have the equation-of-state 
parameter of the cosmological constant. In this case 



E(z) = y/n m (l + z) 3 +Q 



A. 



(7) 



with Sid = ^a- 

The linear growth function, in Eq. is defined 

as D(z) = 5 m (z)/8 m (z = 0) and it is obtained as a 
solution from the following equation (see |Pace, Waizmann 
& Bartelmann||2010 for details) 



3 n n 



2 a 5 E 2 



0, 



(8) 



where the derivatives are taken in relation to the scale factor 
a. 

On the other hand, the equation governing the total gas 
mass (p g ) in the haloes is 



a(M, z) = a{M, z = Q)D(z), 



(2) 



is the linear rms density fluctuation in spheres of comoving 
radius R containing the mass M, and D(z) is the linear 
growth function. 

The density of baryons is proportional to the density 
of dark matter if we consider that the baryon distribution 
traces the dark matter. Thus, the fraction of baryons at 
redshift z that are in structures is given by (see, e.g 
et al.|2006||Pereira fc Miranda|2010 | 



Daigne 



Ps = 



d 2 AL d AL 



dVdt dVdt 



+ a b (t). 



(9) 



The first term on the right hand side of equation |9| 
represents the stars which are formed by the gas contained 
in the haloes. Using a Schmidt law ( Schmidt] 1959, 19631 we 
can write for the star formation rate 



d l M* 
dVdt 



(t) = k Ps (t), 



(10) 
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where k is the inverse of the time-scale for star formation. 
That is, k = l/n. 

The second term on the right hand side of equation (|9| 
considers the mass ejected from stars through winds and 
supernovae. Therefore, this term represents the gas which is 
returned to the 'interstellar medium of the system'. Thus, 
we can write (see, e.g., |Tinsley|1973| ) 



dVdt 



(m — m r )$(m)\t(i — T m )dm, (11) 



m(t) 



where the limit m(t) corresponds to the stellar mass whose 
lifetime is equal to t. In the integrand, m r is the mass of the 



remnant, which depends on the progenitor mass (see Tinsley 
1973| for details), and the star formation rate is taken at the 
retarded time (t — r m ), where r m is the lifetime of a star of 
mass m. 

For all stars formed in the haloes, it is used the 
metallicity-independent fit of |Scalo| ( [T986| ; |Copi| ( |1997[ ) 



M 

M© 



+ 



lop 



M 
M© 



log 10 (r m ) = 10.0-3.6 lo gl 

\ .\ i . / \ ,\ i . / 

(12) 

where r m is the stellar lifetime given in years. 

In equation the term $(m) represents the initial 

mass function (IMF) which gives the distribution function 
of stellar masses. Thus, 



$(m) = Am 



-(l+x) 



(13) 



where x is the slope of the IMF, and A is a normalization 
factor determined by 



m$(m)dm = 1. 



(14) 



O.IMq 



Numerical integration of Q produces the function p g (t) 
at each time t (or redshift z). Once obtained p s (i), we return 



to Eq. ( 10 1 in order to obtain the cosmic star formation rate 



(CSFR). Just replacing *(i) by p\(t), we have 



kp g 



(15) 



It is worth stressing that although we did not take 
into account the stellar feedback processes on the derivation 
(see, 



e.g. 



Christensen et al. 2010 for this 



of the CSFR 

issue), our models, as we will discuss below 



have good 

agreement with observational data taken from [Hopkins 



(2004 20071 at lower redshifts (z < 5). Furthermore, the 



CSFR obtained in the present work has good agreement 
with that one derived by [Springel fc Hernquist| ( |2003[ ) from 
hydrodynamic simulations. Regardless, we should comment 
on this limitation of the model in its present form. In 
particular, stellar feedback processes can modify the time- 
scale (t s ) for star formation, for example, through radiation, 
winds, and supernova events from massive stars. As a main 
result, the star formation efficiency, e*, embedded in the 
normalization of the CSFR at z 
2010 for this issue) and t, 
and e 



good agreement with observations at z < 2, where data are 
less scattered, the values used in this work should represent 
reasonable mean values of these quantities over the whole 
interval [0, 2i n i]- Certainly, the inclusion of stellar feedback 
processes would be an interesting refinement to introduce in 
future works. 



2.2 the dark-energy models and the input 
parameters 

The last point we have to consider for the characterization 
of the CSFR is the dark-energy component of the Universe 
through its equation-of-state parameter w(a). In the present 
work, we consider two cases: cosmological constant where 
w(a) = — 1 and the Chaplygin gas. 

In particular, the Chaplygin gas is characterized by 
a fluid with an equation-of-state p — ~A/p a . This dark 
energy fluid has been tested against observational data 
(e-g, 



as, for example, SNIa 
cosmic microwave background (e.g 
power spectrum (e.g. 
and it configures in a 
cosmological constant. 



Colistete Jr. & Fabris 



Piattella 20101, and 



Fabris, Velten & Zimdahl 20101 



2005), 



strong alternative candidate to the 
Its equation-of-state parameter is given by 



w(a) 



A 

~ A + B a-'^+V 



(16) 



In Eq. (161, the constants A and B are (see, e.g., Pace, 
Waizmann &: Bartelmann|2010 for details) 



A = —w {£ldPc) 1+a and B = (l + w )(fl d p c ) 1+a . (17) 

The present value of the equation-of-state parameter is 
related to A and B by 



A 



wo = 



A + B 



(18) 



The cosmological parameters we have used in this work 
are: Q d = 0.762, fi m = 0.238, D. h = 0.042, and Hubble 
constant H = 100/ikms -1 Mpc~ x with h = 0.734. For the 
variance of the overdensity field smoothed on a scale of size 
we consider as = 0.8. The parameters associated 
to the Chaplygin gas are a — 1.0 (classical Chaplygin gas) 
and a — 0.2 (generalized version). In both cases we consider 
w = -0.8. 

In the column 1 of Table 1 is shown the name of the 
models. In column 2 we present the slope of the initial mass 



function (x in equation 13 1, the time-scale for star formation 
is presented in column 3, the redshift (z p ) where the CSFR 
peaks is presented in column 4, and finally in column 5 we 
have the kind of dark fluid. 

All models presented in Table 1 have good agreement 
with observational data. In particular, it was performed x 2 
analysis over these models, obtaining the reduced chi-square 
defined as Xr = X 2 /dof (where "dof means "degrees of 
freedom"). All of these models satisfy Xi < 1. In Fig. 1 



(see Pereira & Miranda we present the CSFR derived from equation ( 15 1 for three 



would be functions of time. In 
not constants can modify the shape of 



principle, r E 

the CSFR at higher redshifts. On the other hand, due to the 



models of Table 1. The observational points are taken from 
Hopkins| ( [2004] [2007] , 

We can see from the results in Table 1 and Fig. 1 that 
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Table 1. The input parameters used to obtain the Cosmic 
Star Formation Rates. All CSFRs have good agreement with 
observational data. The redshift z- la i associated with the 
beginning of star formation is 20. In the fifth column, A represents 
the cosmological constant and a represents the Chaplygin gas. 



CO 
I 



O 
CL 



© 



CSFR x (IMF) t s Gyr 



dark fluid 



Al 


1.35 


2.0 


3.54 


A 


A2 


1.35 


3.0 


2.94 


A 


A3 


0.35 


1.0 


3.29 


A 


Al 


1.35 


2.0 


2.75 


a = 0.2 


A5 


1.35 


3.0 


2.21 


a = 0.2 


A6 


0.35 


1.0 


2.52 


a = 0.2 


AT 


1.35 


2.0 


2.42 


a = 1.0 


A8 


1.35 


3.0 


1.91 


a = 1.0 


A9 


0.35 


1.0 


2.21 


a = 1.0 



0.1 



0.01 



0.001 




The connection between the past CSFR and the rate of 
binary merger is given by p* c (z) through the relation 



P*c{z) = 



t(z) 



(1 + Zf) 



P(U)dt d 



(20) 



where p*(zf) is the CSFR obtained from Eq. (151, P(td) 
is the probability per unit of time of merging after the 
formation of the progenitor, including both the evolutionary 
time for the formation of the compact binary and the time 
for the compact binary to coalesce, and the (1 + Zf) term 
in the denominator considers the time dilatation due to the 
cosmic expansion. 

The time delay td makes the connection between the 
redshift z at which a compact binary system mergers, and 
the redshift Zf at which its progenitor was formed. As 



discussed by Regimbau Regimbau & Hughes ( 2009 1 , it can 
be calculate by 



1 



dz' 



{l + z')E{z'Y 
The probability P(td) is described of the form 



P(td) oc 



(21) 



(22) 



As mentioned by Regimbau & Hughes (2009), this 
form accounts for the wide range of merger times observed 
in binary pulsars. This form was also used by de Freitas 
Pacheco (1997), [Regimbau fc de Freitas Pacheco ( 2006 1 , and 



define -P(id) as 



de Freitas Pacheco et al. (20061 in their works. Thus, we 



Figure 1. The CSFR derived in this work compared to the 
observational points (HP) taken from |HopkTns] ( |2004| |2007[ l . The 
models are described in Table 1. 



the cosmological constant produces higher amplitudes than 
the Chaplygin gas. That means, the process of baryonic 
matter infall from the haloes is more efficienty, for the same 
set of parameters, if the dark energy fluid is the cosmological 
constant. Another characteristic which can be seen from 
Fig. 1 and Table 1 is that the Chaplygin gas decreases 
the redshift where the CSFR peaks when compared to the 
cosmological constant. 



P{U) = 7- 



(23) 



where B is a normalization constant, and the probability 
function P{td) is normalized in the range of to — 15 Gyr for 
some minimal delay To- Therefore, 



15Gyr 



B 

U 



dtd = 1. 



(24) 



Specifically, we consider that to = 20Myr for NS-NS 
systems, to = lOMyr for NS-BH systems, and to = 100 Myr 
for BH-BH systems ( |Bulik, Belczynski fc Rudak|2004 l. 

With these assumptions, the merger rate per unit of 
redshift can be written as 



3 THE COALESCENCE RATES 

We assume that the coalescence rates track the CSFR but 
with a delay id between the formation of the binary system 



and the final merger (see Regimbau fc Hughes|2009[ ). Thus, 
we can write 



p°(z) = p°(0) x 



P*c(Z) 



P*c 



(0)' 



(19) 



where Pc(z) is the rate at which binary systems are observed 
to merger at redshift z, and p°(0) is the same rate in our 
local Universe. 



dR° c .o dV 
-dz- =P ^ Z) lz-> 



(25) 



where dV is the comoving volume element which is given by 

.2 c dz 



dV = 4Tvr(zY 



H E{z) 



dz 
Ejzj- 



(26) 



In Eq. (261, r(z) is the proper distance, whose 
expression is 



(27) 
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Figure 2. The cosmic coalescence rates for NS-NS binaries. The 
models consider to = 20 Myr. These models are obtained from 
the CSFR presented in Fig. 1 (see also Table 1 for the main 
parameters). 



Figure 4. The cosmic coalescence rates for BH-BH binaries. The 
models consider to = 100 Myr. These models are obtained from 
the CSFR presented in Fig. 1 (see also Table 1 for the main 
parameters) . 



o 



3.5 
3 
2.5 



•a. 1 



1-5 



o 
•a. 



1 

0.5 






20 



Figure 3. The cosmic coalescence rates for NS-BH binaries. The 
models consider to = 10 Myr. These models are obtained from 
the CSFR presented in Fig. 1 (see also Table 1 for the main 
parameters). 



Note that the expansion function E(z) is dependent on 
the kind of dark-energy fluid as shown by Eq. (fsj) . Then, 
using the formalism described in this Section, it can be 
determined the cosmic coalescence rates for NS-NS binaries, 
NS-BH, and BH-BH systems up to redshift z ~ 20. Figure 
2 presents the cosmic coalescence rates, normalized to the 
local value p°(0), for NS-NS binaries. On the other hand, 
Figures 3 and 4 respectively show the cosmic coalescence 
rates for NS-BH and BH-BH binaries. 

As expected, due to the behavior of the CSFR, 
the amplitudes of the coalescence rates produced by the 
Chaplygin gas cosmology are lower than those produced 
by the cosmological-constant cosmology. In particular, for 
the cosmological constant (model Al) p°(z)/p°(0) reaches a 
maximum amplitude at redshift z = 2.27 for NS-NS binaries, 
at z = 2.45 for NS-BH and at z = 1.86 for BH-BH systems. 
In the case generalized Chaplygin gas (a — 0.2 — model A4) 



we note that p°(z)/p°(0) peaks at z = 1.63 for NS-NS, at 
z = 1.81 for NS-BH and at z = 1.31 for BH-BH. The last 
case, classical Chaplygin gas (a = 1.0), reaches a maximum 
amplitude at z = 1.34 for both NS-NS and NS-BH systems. 
On the other hand, for BH-BH binaries the maximum value 
of the coalescence rate is reached at z = 1.13. Thus, the 
position of the peak of the coalescence rate is dictated by 
the value of t and also by the kind of dark-energy fluid. 

As we will see in the next Section, the different 
behaviors for the coalescence rates produced by different 
dark-energy cosmologies will produce different values for 
the signal-to-noise ratios of advanced LIGO and Einstein 
Telescope. 



4 GRAVITATIONAL WAVE BACKGROUND 
4.1 Compact binary systems 

The spectrum of a stochastic background of GWs is 
characterized by the closure energy density per logarithmic 



frequency spam, which is given by (see, e.g., Allen 1997 
| Allen fc Romano|1999[ ) 



dpG\ 



, , (28) 
p c d log v ohs 

where pgw is the gravitational energy density, and v h s is 
the frequency in the observer frame. 



The above equation can be written as (see, e.g., Ferrari, 
Matarrese fc Schneider|1999 1 



C^pc 



(29) 



where F„ obB is the gravitational wave flux (given in 
erg cm" 2 Hz" 1 s" 1 ) at the observer frequency v a bs integrated 
over all cosmological sources. Therefore, 



f, obs dR° c (z). (30) 
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See that dR®(z)/dz is the merger rate per unit of 
redshift (Eq. 251. In order to solve Eq. (301, it is needed 
to determine the gravitational wave fluence {fv obs ), in 
the observer frame, produced by a given compact binary 
coalescence. Follo wing |Regimbau fe de Freitas Pacheco 



(2006) 



Regimbau (20111); Zhu et al. (2011b; IMarassi et al 



( 2011[); |Rosado| j2011| >; |Wu, Mandic fc Regimbau| ( |2012| ); 
Regimbau (20121, f VohB can be written as 



1 dE G \ 



dv 



4ndl 



(31) 



where c!l = r ( 2 )(l + z ) is the luminosity distance, r(z) is 
the proper distance (see Eq. 27 1, dEGw/du is the spectral 
energy, and v = ^ b s (l + z) is the frequency in the source 
frame. 

In the quadrupolar approximation, the spectral energy 
emitted by a compact binary system, with masses mi and 



rri2, which inspirals in a circular orbit is given by (Peters & 



Mathews 1963) 



a 



3.7e-09 
3.0e-09 
2.2e-09 
1 .5e-09 
7.5e-10 
0.0e+00 



Al A4 


- A7 

- 


/ /<" — ->0\ 




/ "" : --\ 




~ // \ 




/ * 




K NS-NS (x = 20 Myr) 
i.i.i 


1 . 



300 



600 900 

v bs ( Hz ) 



1200 



1500 



Figure 5. Spectrum of the gravitational energy density 
parameter Hgw- The results are shown for double neutron stars 
(NS-NS). 



dEp\ 
dv 



Kv 



-1/3 



where 



K 



(Gtt 



,2/3 



3 (mi + m2) 1 ' 3 



(32) 



(33) 



It will be considered that the GW background has the 
value of the maximum frequency limited by the 'last stable 



orbit' (LSO). Then, following Sathyaprakash (20011 



= = 1.5 



M 



2.8M Q 



kHz, 



(34) 



where M is the total mass of the system (M = mi + ma). 

In the present study, we consider mi = m.2 = 1.4 Mq 
for NS-NS binaries, while for NS-BH are used mi = 1.4 M Q 
and wi2 = 7.0 Mq. Concerning to BH-BH systems we have 
used mi = m2 = 7.0 Mq. With these considerations, the 
maximum frequency is ^lso = 1.5 kHz (500 Hz) for NS-NS 
(NS-BH). For BH-BH binaries we have u LSO = 300 Hz. 

There is one last point to consider before calculating the 
spectrum of the stochastic background of GWs. This point is 
related to the value of the local merger rate per un it volume 
Pc(0). As discussed by Regimbau & Hughes ( |2009 1, the local 
merger is usually extrapoled by multiplying the rate in the 
Milk Way with the density of equivalent galaxies. 

Current estimates give p°(0) = (0.01 — 10) x 
Myr^Mpc" 3 for NS-NS, an d pg(0) = (0.001 - 1) x 



2009 



Myr^Mpc" 3 for NS-BH (see |Regimbau fc Hughes 
and the references therein). In the present work, it is 
considered p°(0) = 1.0 Myr^Mpc" 3 for NS-NS, p°(0) = 
O.lMyr^Mpc" 3 for NS-BH, andp°(0) = 0.01 Myr^Mpc -3 
for BH-BH. 

Thus, using the formalism above, we can obtain the 
characterization of the stochastic background of GWs 
formed by the coalescence of compact binary systems. In 
particular, the Figures 5, 6, and 7 present the spectra of 
the gravitational energy density parameter JIgw versus the 
observed frequency f bs. 

The density parameter increases as !/^ 3 at low 



a 



7.5e-10 



6.0e-10 - 



4.5e-10 



3.0e-10 



1.5e-10 



0.0e+00 




200 300 

v bs ( Hz ) 



500 



Figure 6. Spectrum of the gravitational energy density 
parameter Oq^t. The results are shown for the coalescence of 
neutron star-black hole (NS-BH) systems. 



frequencies and reaches a maximum amplitude ~ 3.8 x 
10 -9 around 375 Hz for NS-NS systems if is used the 
CSFR-A1 (cosmological constant). It is worth stressing 
that calculations performed by [Regimbau fc de Freitas] 
Pacheco] ([2006), using Monte Carlo methods for obtaining 
the coalescence rates, produced similar results. 

In particular, these authors obtained maximum 
amplitude of about 1.1 x 10" 9 around 670 Hz for NS- 
NS binaries (considering their fiducial CSFR). However, 
observing the distribution of coalescences as a function of 
the redshift derived by Regimba u fc de Freitas Pacheco] 



(20061 we note that it peaks at z ~ 1.5. On the other hand, 
our coalescence rates peak at z ~ 2.9 — 3.5 (models Al to 
A3 which correspond to the cosmological constant as dark- 
fluid). In this way, the maximum value of f^Gw is shifted 
to lower frequency than that obtained by [Regimbau fc de| 
|Freitas Pacheco| ( |2006| ). 

In order to assess the detectability of a gravitational 
wave signal, one must evaluate the signal-to-noise ratio 
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Figure 7. Spectrum of the gravitational energy density 
parameter f^GW- The results are shown for the coalescence of 
black hole-black hole (BH-BH) systems. 



(S/N), which for a pair of interferometers is given by (see, 
for example, |Christensen||1992| |Flanagan||l993"l |Allen||1997l 



de Araujo, Miranda fc A guiar 2002] |2004| |Regimbau fc de| 
Frcitas Pacheco||2006 1 



(S/N) 



50tH 



T 



dv- 



(2), 



(35) 



where S^ is the spectral noise density, T is the integration 
time, and j(u) is the overlap reduction function, which 
depends on the relative positions and orientations of the two 
interferometers. For the 7(f) function we refer the reader to 
Flanagan ( 1993 1 who was the first to calculate a closed form 



for the LIGO observatories. 

The noise power spectral density of the advanced LIGO 
can be found at LIGO website. However, we used here the 



analytical fit given by Mishra et al. (20101. Its expression is 



Sh(u) = So 



10 



16-40-7. 9) 2 



+ 2.4 x 10" 



- 50 + 0.08x- 4 '' 



+ 123.35 

= 00 if v < v, 



0.23 x 2 + 0.0764 a; 4 
1 + 0.17a; 2 



if v > v s 



(36) 



where x = v/vq, v stands for the frequency, vq = 215 Hz, 
So = Hn 49 Hz-\ and u s is a low-frequency cut-off that can 
be varied, and below which Sh(f) can be considered infinite 
for all practical purposes (here, we choose v s = 10 Hz). 

On the other hand, a possibility for a third generation 
ground-based gravitational wave detector is the Einstein 
Telescope (ET). The basic design of this inteferometer is 



curvetp 




Hild, Chelkowski & Freise 


2008 


Hild et al. 


2010 


Punturo et al.|2010a|b Sathyaprakash et al.|2012i. 



1 See alsolhttp: / /www.et-gw.eu/etsensitivitiesl In particular, the 
ET-C and ET-D sensitivity curves correspond to a xylophone 
configuration consisting of a pair of detectors. The first one 
operating at low-frequency (1 — 100 Hz) and the second detector 
operating at high-frequency (100— ~ 10 kHz). 



Table 2. The main characteristics of the models NS-NS and their 
respectives signal-to-noise ratios (S/N) for a pair of 'advanced 
LIGOs' and Einstein Telescope in triangular configuration. The 
integration time is T = 1 yr. We also show the redshift zrjc at 
which the duty cycle becomes equal to 1 (transition between the 
popcorn and the continuous stochastic regime). 



NS-NS 


^GWmax 


v P (Hz) 


(S/N) 


(S/N) 


2DC 








Adv. LIGO 


ET 


D = 1 


Al 


3.84 x 10~ 9 


375 


1.46 


335 


0.55 


A2 


3.11 x 10~ 9 


415 


1.11 


256 


0.57 


A3 


3.57 x 10~ 9 


390 


1.32 


305 


0.56 


A 1 


2.70 x 10~ 9 


446 


0.92 


211 


0.61 


A5 


2.27 x 10~ 9 


485 


0.71 


168 


0.64 


A6 


2.53 x 10~ 9 


461 


0.84 


194 


0.62 


A7 


2.34 x 10~ 9 


482 


0.76 


174 


0.64 


A8 


2.00 x 10~ 9 


521 


0.61 


142 


0.67 


A9 


2.21 x 10~ 9 


498 


0.70 


161 


0.65 



Here, we use the ET-B sensitivity curve of |Hild,| 



Chelkowski & Freise ( 2008 1 with an analytical fit taken from 
Mishra et al.| < |2010[ >. 



= So [ 



fol , bo 1 b^ , b 

a\x + a^x + asx + a^x 



if v < v s 



vjvo, v stands for the frequency, vo 
-1 



if v v s 

(37) 
= 100 Hz, 



where x 

So — 10 -50 Hz _i , and v s is a low-frequency cut-off that can 
be varied, and below which Shiy) can be considered infinite 
for all practical purposes (here, we choose v s — 10 Hz). The 
coefficients in Eq. ([371 have the values 



a 2 
a 3 

0,4 



2.39 x 10" 
0.349, 
1.76, 
0.409, 



bi = 


-15.64, 


62 = 


-2.145, 


63 = 


-0.12, 


b 4 = 


1.10. 



(38) 



We consider that the ET has a triangular configuration 
( Hild et al.||2010 1 with an overlap reduction function given 
by Regimbau (20111. In Tables 2, 3, and 4 we summarize 
the main characteristics of the models. We show the values 
of the maximum amplitude (^Gw max ) of the stochastic 
background, the frequency (u p ) where Qgw peaks, and 
the signal-to-noise ratios for advanced LIGO and ET 
interferometers. Note that for all kind of compact binaries 
we have (S/N) < 1 for a pair of advanced LIGOs. On the 
other hand, for ET in triangular configuration should be 
possible, in principle, to obtain high values for (S/N). 

Concerning the nature of the GW background, it is 
determined by the duty cycle which is defined as the ratio, in 
the observer frame, of the typical duration of a single burst 
t to the average time interval between sucessive events (see 
Regi mbau fc de Freit as Pachcco 2006; Regi mbau fc Hughes| 
20091 



(39) 



where 
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Table 3. The main characteristics of the models NS-BH and their 
respectives signal-to- noise ratios (S/N) for a pair of 'advanced 
LIGOs' and Einstein Telescope in triangular configuration. The 
integration time is T = 1 yr. We also show the redshift zrjc at 
which the duty cycle becomes equal to 0.1 (transition between 
the shot noise and the popcorn regime). 



NS-BH 


^GWmax 


v P (Hz) 


(S/N) 


(S/N) 


ZDC 








Adv. LIGO 


ET 


D = 0.1 


Al 


6.53 X 10~ 10 


123 


0.52 


116 


0.93 


A2 


5.25 x 10~ 10 


136 


0.40 


88 


0.99 


A3 


6.06 x 10~ 10 


128 


0.47 


105 


0.94 


A 1 


4.53 x 10~ 10 


146 


0.33 


73 


1.15 


A5 


3.79 x 10~ 10 


159 


0.26 


58 


1.28 


A6 


4.25 x 10~ 10 


151 


0.30 


67 


1.19 


A7 


3.92 x 10~ 10 


158 


0.27 


61 


1.30 


A8 


3.35 x 10~ 10 


172 


0.22 


49 


1.51 


A9 


3.69 x 10~ 10 


164 


0.25 


56 


1.37 



Table 4. The main characteristics of the models BH-BH and their 
respectives signal-to- noise ratios (S/N) for a pair of 'advanced 
LIGOs' and Einstein Telescope in triangular configuration. The 
integration time is T = 1 yr. 



BH-BH 


^GWmax 


v v (Hz) 


(S/N) 


(S/N) 








Adv. LIGO 


ET 


Al 


1.68 x 10" 10 


81 


0.17 


36 


A2 


1.40 x lO" 10 


89 


0.14 


29 


A3 


1.58 x 10~ 10 


81 


0.16 


33 


Al 


1.21 x 10~ 10 


96 


0.11 


24 


A5 


1.04 x lO" 10 


103 


0.09 


20 


A6 


1.14 x lO" 10 


99 


0.10 


22 


A7 


1.06 x 10~ 10 


103 


0.09 


20 


A8 


9.21 x 10~ n 


110 


0.08 


17 


A9 


1.00 x lO" 10 


106 


0.09 


19 



5c 5 



256tt 8 / 3 G 5 / 3 



( 1 + z ')m c ]- 5/3 K 8/3 , 



(40) 



with /l being the lower frequency bound of the detector, 
and m c represents the chirp mass which is given by 



(mim 2 ) 3/5 



(mi + m 2 



,1/5 



(41) 



Figure 8 presents D(z) for the NS-NS binaries while 
Fig. 9 shows the duty cycle for NS-BH systems. In Figure 
10 we have the duty cycle for BH-BH systems. In these plots 
we have considered /l = 10 Hz. 

Concerning the duty cycle, there are three different 
regimes for this parameter (see, e.g., Rosado 2011 Wu, 



Mandic fc Regimbau|[2012| |Regimb au 2012). The first case 



appears when D(z) < 0.1. In this case, we have the so-called 
'shot noise regime' consisting of a sequency of widely spaced 
events. That means the sources can be resolved individually. 

The second case appears when 0.1 < D(z) < 1. We have 
the 'popcorn noise regime' That means the time interval 
between two sucessive events could be closer to the duration 




Figure 8. The duty cycle as a function of the redshift for NS-NS 
binaries. We consider /l = 10 Hz. The horizontal line at D(z) = 1 
represents the transition between the popcorn regime and the 
continuous stochastic background. 



of a single event. In reality, near to D(z) = 1 the events may 
overlap, making it difficult to identify individual events. 

The third case appears when D(z) 1. In this case, 
we have a 'continuous background'. The signals overlap to 
produce a continuous stochastic background. 

In Table 2 we include the value of the redshift zuc at 
which the background becomes continuous (D(z) > 1). In 
Table 3 we present the redshift of transition between the shot 
noise and the popcorn regime (D(z) > 0.1). In particular, 
for NS-NS binaries a continuous background is established 
for sources situated at cosmological distances z ~ 0.5 — 0.6 
(z ~ 0.6 — 0.7) for the cosmological constant (Chaplygin gas) 
cosmology. 

On the other hand, for the NS-BH binaries the signals 
change from shot noise to popcorn regime at z ~ 0.9 — 1.0 
(z ~ 1.0 — 1.5 for Chaplygin gas). For the BH-BH systems 
the signals are always within the shot noise regime. However, 
note that reducing the values of the local merger p®(0) in 
relation to those values used in the present work, it would 
reduce the values of the (S/N) besides changing the regimes 
(or values of D(z)) of the gravitational wave backgrounds. 
In particular, see that (S/N) and D(z) are proportional to 

p2(o). 



4.2 Collapse of stars to form black holes 

In order to determine the background of GWs generated by 



stars which collapse to black holes, we re- write Eq. ( 30 1 as 



F^s = 
where now we have 

di?BH 

dz 

and for fu ohB we have 



Jo 



cIRbh(z), 



= pi,(z)$(m)dj 



dV_ 

dz ' 



(42) 



(43) 
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core collapse to form BH 



100 

v bs ( Hz ) 



1000 



Figure 9. The duty cycle as a function of the redshift for NS- 
BH binaries. We consider /l = 10 Hz. The horizontal line at 
D(z) = 0.1 represents the transition between the shot noise and 
the popcorn regime. 



Figure 11. Spectrum of the gravitational energy density 
parameter Qqw We consider an efficiency of generation of 
gravitational waves £gw = 10~ 4 . 




Figure 10. The duty cycle as a function of the redshift for BH- 
BH binaries. We consider /l = 10 Hz. These compact binaries 
produce a signal of the kind shot noise. 



7TC 2 



(44) 



1987) 



The dimensionless amplitude /ibh is given by ( |Thorne| 



ft BH * 7.4 x 10- 20 4^ 4 M _ 



© 



\lMpc 



(45) 



where £gw is the efficiency of generation of GWs, and m r is 
the mass of the black hole formed. 



It is worth mentioning that Eq. ( 45 1 refers to the black 



hole 'ringing', which has to do with the de-excitation of the 
black hole quasi-normal modes. 

The collapse of a star to black hole produces a signal 
with frequency i/ Q bs given by 



^obs 



1.3 x 10 4 Hz 



MfT 







(46) 



We will consider that black holes are formed from stars 
with 40 Mq ^ m 5j 140 Mq. The lower limit is consistent 
with recent results derived from the X-ray pulsar CXO 
J164710. 2— 455216 which show that the progenitor to this 
pulsar had an initial mass ~ 40 M ( |Muno et al.|2006[ ). On 
the other hand, the mass of the black hole remnant is taken 



to be the mass of the helium core before collapse (see Heger 
fc Woosley|2002) . Thus, 



13 

m r = m Hc = 7^( m - 20 M ©)- 



(47) 



With these considerations, we can obtain the spectrum 
of GWs produced by cosmological black holes. The Fig. 
11 shows the spectrum of the gravitational energy density 
parameter Ogw as a function of the observed frequency v b s 
for the models with the highe st (S/N) of Table 5. These 
curves consider egw = 10 -4 (Loftier, Rezzolla & Ansorg 



2006) 



We can see that the spectra peak at Hgw ~ 3 x 10 — 
10~ 7 dependent on both CSFR parameters and dark-energy 
component. Note that only two models have transition from 
shot noise to popcorn regime (which corresponds to D — 
0.1). All the other models correspond to shot noise signals. 

The Fig. 12 presents the duty cycle generated by the 
collapse of stars to form black holes. Note that in this case 
we calculate the duty cycle as 



D(z)= / f(l + z 



dRsm 
dz< 



(48) 



with f = 1ms (Ferrari, Matarrcs e fc Schneider|1999[ ). 

We can observe from Table 5 that the frequency 
where the spectra peak neither depend on the dark- 
energy cosmology nor on the CSFR parameters. Another 
characteristic of this kind of source is the very high (S/N) 
produced for ET. 
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Table 5. The main characteristics of the models 'core-collapse to 
form black holes'. The efficiency of generation of GWs is £qw 

1 The signal-to -noise ratios (S/N) for a pair of 'advanced 
LIGO' are determined for an integration time T = lyr. For ET 
we consider triangular configuration. We also show the redshift 
2DC a t which the duty cycle becomes equal to 0.1 (transition from 
shot noise to popcorn regime). 



BH 


^GWmax 


v p (Hz) 


(S/N) 


(S/N) 


ZDC 








Adv. LIGO 


ET 


D = 0.1 


Al 


6.80 X 10~ 9 


200 


0.14 


152 




A2 


5.00 x 10~ 9 


200 


0.08 


96 




A3 


9.02 x 10" 8 


200 


1.76 


1900 


3.38 


Al 


4.70 x 10" 9 


200 


0.07 


86 




A5 


3.63 x 10" 9 


200 


0.04 


58 




A6 


6.29 x 10" 8 


200 


0.85 


1100 


6.92 


A7 


4.06 x HT 9 


200 


0.05 


68 




A8 


3.20 x 10" 9 


200 


0.03 


48 




A9 


5.46 x 10" 8 


200 


0.62 


870 






Figure 12. The duty cycle as a function of the redshift for the 
cosmological Population of black holes. The horizontal line at 
-D(z) = 0.1 represents the transition from shot noise to popcorn 
regime. 



It is worth stressing that Zhu, Howell & Blair ( 2010 1 has 



recently estimated the upper limit on the energy density, 
of a stochastic gravitational wave background, produced 
by the core-collapse supernovae leaving black holes as 
remnants. The authors showed that considering Gaussian 
source spectra would be possible to detect GW signals 
with £ G w ~ 10~ 5 (10~ 7 ) for Adv. Ligo (ET). Another 
work centerded on the GW backgrounds produced by core- 
collapse supernovae of Population III and Population II 



stars was developed by Marassi, Schneider & Ferrari (2009). 
In particular, the authors studied the cosmic transition of 
Population III to Population II using waveforms derived 
from recent 2D numerical simulations. The GW efficiencies 
used by these authors were £gw ~ 10~ 7 for Population II 
progenitors and egw ~ 10 for Population III progenitors 
(with initial masses ranging between 100 Mq to 500 Mq). 
Here, we have adopted the GW spectrum of |Thorne| {1987), 
with efficiency egw ~ 10 -4 because the core-collapse energy 



spectrum will affect the results for both cases, A-CDM and 
Chaplygin gas, exactly the same way. 



5 COSMOLOGICAL SPECTRUM PRODUCED 
BY ALL SOURCES 

It is worth stressing that a stochastic background of 
gravitational waves is expected to arise from a superposition 
of a large number of gravitational wave sources of 
astrophysical and cosmological origin. In particular, the 
Wilkinson Microwave Anisotropy Probe (WMAP) results 
suggest an early epoch for the reionization of the Universe 
(see, e.g., Jarosik et al. 20111. In this way, a pre-galactic 



Population should be formed at high redshifts to account 
for these results. The cosmic star formation history is 
determined by the interplay between the incorporation of 
baryons into collapsed objects and return of baryons into 
diffuse state (e.g., gaseous clouds). 

Thus, the formation of different objects as, for example, 
NS-NS binaries, NS-BH binaries, BH-BH systems, core- 
collapse to form black holes, among others, are directly 
related to the CSFR. On the other hand, different dark- 
energy scenarios could give different signatures on the 
background through the expansion function E{z) (Eq.jfjj). As 
the gravitational wave background could trace the behavior 
of the Universe up to redshift ~ 20, then it should be 
possible to infer if there is a temporal dependence of the 
dark-energy equation of state. That is, if w(a) 7^ 0. 

In this way, the detection and characterization of a 
stochastic background of GWs could be used as a tool for 
the study of the Universe at high redshifts. In particular, the 
gravitational wave signals produced at different cosmological 
distances by the sources discussed above could overlap at 
a given frequency j/ b s to produce a stochastic background 
over a large range in frequency. 

In Figure 13 we show the collective contribution of the 
three compact binary sources investigated here. In Tabic 
6 we summarize the main characteristics of these models. 
Although the (S/N) of the collective spectra are dominated 
by the NS-NS binaries, we note that the NS-BH and BH- 
BH binaries pull the peak of the collective spectra for lower 
frequencies than those observed if we only consider the NS- 
NS systems. 

In Figure 14 we include the core-collapse supernovae 
together the compact binary systems in the calculation of 
the gravitational wave spectra. We can see that same for a 
low efficiency of generation of gravitational waves (egw = 
10 -4 ) the collective spectra show a clear signature of this 
kind of source when compared to the spectra only derived 
with compact binaries. 

In Table 7 we present the main caractheristics of the 
collective spectra with all sources studied in this work. 

We note that core-collapse supernovae have an 
important contribution for the shape of the collective 
spectrum for frequencies 60— ~ 300 Hz. In particular, the 
frequency where the collective spectra peak is completely 
dominated by the core-collapse supernovae. See that model 
A3 could produce (S/N) ~ 3 (~ 2200) for advanced LIGO 
(ET) in the cosmological-constant cosmology. On the other 
hand, model A6 could produce (S/N) ~ 2 (~ 1300) for 
advanced LIGO (ET) in the Chaplygin gas cosmology. In 
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100 
v obs ( Hz ) 

Figure 13. Collective spectra of the three compact binary sources 
studied in this work. The curves represent the models Al, A4, and 
A7. 



Table 6. The main characteristics of the collective contribution 
of the three compact binary sources, namely, NS-NS binaries, NS- 
BH systems and BH-BH binaries. The signal-to- noise ratios (S/N) 
for a pair of 'advanced LIGO' are determined for an integration 
time T = lyr. We also present the (S/N) for ET in triangular 
configuration. 



Model 


^GWmax 


v P (Hz) 


(S/N) 


(S/N) 








Adv. LIGO 


ET 


Al 


4.06 X 10~ 9 


321 


2.15 


486 


A2 


3.27 x 10~ 9 


346 


1.64 


373 


A3 


3.77 x 10" 9 


331 


1.95 


442 


Al 


2.81 x 10" 9 


374 


1.35 


308 


A5 


2.34 x 10~ 9 


398 


1.08 


246 


A6 


2.63 x 10~ 9 


385 


1.24 


283 


AT 


2.42 x 10" 9 


396 


1.12 


254 


A8 


2.04 x 10" 9 


413 


0.91 


207 


A9 


2.27 x 10" 9 


404 


1.03 


235 




10" 



10" 



all cosmological sources 



10 



100 

•'obs 



1000 



(Hz) 



Figure 14. Spectrum of all sources studied in this work. 
Concerning to the core-collapse to form black holes, we use 
e GW = 10 — 4 . The curves represent the models Al, A4, and A7. 



Table 7. The main characteristics of the collective contribution of 
all sources studied in this work. The signal-to- noise ratios (S/N) 
for a pair of 'advanced LIGO' are determined for an integration 
time T = lyr. We also present the (S/N) for ET in triangular 
configuration. 



Model 


^GWmax 


k p (Hz) 


(S/N) 


(S/N) 








Adv. LIGO 


ET 


Al 


1.06 X 10" 8 


200 


2.23 


590 


A2 


8.00 x 10~ 9 


200 


1.69 


433 


A3 


9.36 x 10- 8 


200 


3.26 


2182 


Al 


7.23 x 10" 9 


200 


1.39 


362 


A5 


5.69 x 10" 9 


200 


1.10 


280 


A6 


6.52 x 10~ 8 


200 


1.81 


1262 


A7 


6.19 x 10~ 9 


200 


1.14 


295 


A8 


4.96 x 10" 9 


200 


0.92 


233 


A9 


5.65 x 10" 8 


200 


1.42 


1005 



principle, detecting stochastic backgrounds of gravitational 
waves with high (S/N) could be possible to infer, for 
example, the behavior of the CSFR at high-redshifts as 
well as if there is a temporal dependence of the dark-energy 
equation of state. 



6 COMPACT BINARIES AND 

CORE-COLLAPSE PARAMETERS: 
INFLUENCE ON THE RESULTS 

In the previous Sections, we have analyzed the main 
characteristics of the stochastic backgrounds produced by 
compact binary systems, core-collapse to form black holes, 
and the composite signal of these cosmological objects. We 
verify that higher signal-to-noise ratios can be produced if 
we consider Einstein Telescope in triangular configuration. 
Although these cosmological sources are connected by the 
CSFR, we know that there are uncertainties in the minimum 
coalescence time-scales of NS-NS, NS-BH, and BH-BH 
binaries. On the other hand, the local coalescence rates of 
these systems can vary up to three orders of magnitude. 
In addition, the minimum mass able to form a black hole 
may vary from ~ 25M© to ~ 40Mq. Thus, in this Section, 
we present an analysis of these uncertainties and their 
influence on the stochastic backgrounds discussed here. We 
also discuss if there is a clear difference between the ACDM 
and the Chaplygin gas which would permit constrain both 
the constant feature (or not) of the dark-energy equation of 
state and the CSFR derived for ACDM and Chaplygin gas. 
In particular, we have analyzed: 

a) The local coalescence rate: This parameter acts like 
an offset and it does not modify the shapes of the spectra. 
Note, however, that (S/N) oc Pc(0) and so our values for 
the signal-to-noise ratios can vary from 0.1 to 10 of those 
listed in Tables 2 — 4. Thus, this parameter can only modify 
the values of the (S/N) as (S/N) = p°(0)/p°(0) u X (S/N) u 
where the subscript u means the values used and derived 
in this work. Note that same in the worst case (Pc(0) = 
0.1 x p°(0) u ), it would be possible to have (S/N) > 10 
(ET) for the composite signals of these three binary sources. 
Concerning to the duty cycle, observe that D(z) is also 
proportional to Pc(0). Thus, the redshifts of transition from 
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Figure 15. Collective spectra of the three compact binaries 
taking into account the uncertainties in the parameters. The black 
area describes all the possible GW signals for the ACDM case. 
The gray area represents the GW backgrounds for Chaplygin gas 
with a = 0.2. 



Figure 16. Collective spectra of all sources studied in this work 
and taking into account the uncertainties in the parameters. The 
black area describes all the possible GW signals for the ACDM 
case. The gray area represents the GW backgrounds for Chaplygin 
gas with a = 0.2. 



popcorn to continuous regimes (and from shot noise to 
popcorn regimes) can change according to the values of 
pg(0); 

b) The minimum stellar mass to form a black hole: This 
parameter basically changes the maximum frequency of the 
background formed by core-collapse. In the Case TTlmm — 
25M©, we obtain = 4.8kHz while for m m i n = 40Mq we 
have i/ max = 1.2kHz. In terms of (S/N), if we change m m i n 
from 40Mq to 25M© the signal-to-noise ratios increase in 5% 
in relation to those values present in Table 5. Looking for 
the results in Table 7, collective contribution of all sources, 
(S/N) increases by 2% (50%) for Adv. LIGO (ET). The 
frequency where Ogw peaks is weakly dependent on this 
parameter in both cases ACDM and Chaplygin gas; 

c) Efficiency of generation of gravitational waves (egw): 
Note that SIgw <x £gw- Thus, same with an efficiency of 
generation of GWs ~ 10~ 5 — 10~ 6 could be possible to have 
(S/N) > 10 for ET in triangular configuration (see Table 5); 

d) Coalescence time-scale of NS-NS: We change this 
parameter from 20 Myr to 100 Myr. As a consequence, the 
coalescence rate peaks at z ~ 1.9(1.30) instead of z ~ 
2.27(1.63) for the ACDM (Chaplygin gas with a = 0.2) 
while the (S/N) of the Table 2 typically decreases by 15%. 
Looking for the collective contribution of all sources in Table 
7 we note that (S/N) decreases by 8% (3%) for Adv. LIGO 
(ET). There is just a slight modification of the frequency 
where flaw peaks. 

e) Coalescence time-scale of NS-BH: We change this 
parameter from 10 Myr to 50 Myr. As a consequence, the 
coalescence rate peaks at z ~ 2.1 (1.50) instead of z ~ 
2.45(1.81) for the ACDM (Chaplygin gas with a = 0.2) 
while the (S/N) of the Table 3 typically decreases by 10%. 
Looking for the collective contribution of all sources in Table 
7 we note that (S/N) decreases by 1% for both Adv. LIGO 
and ET. 

f) Coalescence time-scale of BH-BH: We change this 
parameter from 100 Myr to 500 Myr. As a consequence, 
the coalescence rate peaks at z ~ 1.25 (0.83) instead of 
z ~ 1.86 (1.31) for the ACDM (Chaplygin gas with a = 0.2) 



while the (S/N) of the Table 4 typically decreases by 20%. 
Looking for the collective contribution of all sources in Table 
7 we note that (S/N) decreases by 0.8% for both Adv. LIGO 
and ET. 

A question could arise about the uncertainties described 
above: Is it possible to have a clear separation of the two 
backgrounds (ACDM and Chaplygin gas cosmologies) or 
the uncertainties listed above produce a superposition of 
these backgrounds? A second question could also arise: Can 
different dark energy scenarios produce distinct signatures 
on the CSFR? In order to answer these questions, we present 
in Figures 15 and 16 the gravitational wave backgrounds 
with the uncertainties discussed above and for the models 
Al and A3 of Table 1. In these Figures, we just keep fixed 
two parameters: egw = 10 -4 and m m i n = 40 Mq. Below 
^obs ~ 1 kHz there is no superposition between the GW 
signals in the case ACDM and Chaplygin gas (a = 0.2), 
same with all the uncertainties in the parameters. However, 
note that the case Chaplygin gas with a = 1 can not be 
separated from a = 0.2. There is a superposition between 
these two Chaplygin models if we take into account all the 
uncertainties discussed above. 

The second point is related to the CSFR. Looking for 
Figure 17, we see that there is no overlap between the 
case ACDM and Chaplygin gas at z > 2 same with all 
the uncertainties in the parameters. The areas defined by 
ACDM and Chaplygin gas cosmology do not overlap in the 
redshift range [2 — 20]. In principle, with observational data 
less scattered in the range z ~ 2 — 5, it would be possible 
to have a better indication of the dark-energy equation of 
state from the observed CSFR. Otherwise, being detected 
a stochastic background of GWs with high (S/N), as in 
the case of ET, we could work with the inverse problem 
reconstructing the CSFR from the observed background. In 
this way, ET could contribute with a better comprehension 
of how star formation is regulated at high redshifts. 
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Figure 17. Possible CSFRs taking into account all the viable 
models studied in this work. The black are represents the family 
of CSFRs for the ACDM cosmology (models Al to A3) while 
the gray area shows the family of CSFRs for the Chaplygin gas 
(models A4 to A6) as dark-energy component of the Universe. 
Note that there is no overlap at z > 2 between these two dark- 
fluids same with all the uncertainties in the parameters. 



7 FINAL REMARKS 

In this work, we have first studied the main characteristics 
of the gravitational wave signals produced by coalescences of 
NS-NS, NS-BH, and BH-BH binaries up to redshift z ~ 20. 
The coalescence rates are obtained from the hierarchical 
formation scenario recently studied by Pcreira fc Miranda| 
( 2010 1. 

In this formalism, the 'cosmic star formation rate' 



(CSFR) is derived in a self-consistent way, considering the 
baryon accretion rate as an infall term which supplies 
the gaseous reservoir in the haloes. However, here we 
modify their model in order to incorporate different dark- 
energy fluids. In particular, we show results considering 
two different dark-energy components of the Universe: the 
cosmological constant and the Chaplygin gas. 

For NS-NS systems, the shape of the spectrum of the 
gravitational energy density parameter (^gw) has good 
agreement with the case studied by [Regimbau fc de Freitas] 
|Pacheco] |2006) who used numerical simulations based on 
Monte Carlo methods. In particular, we have obtained 
signal-to-noise ratios (S/N) ~ 1.5 for NS-NS, (S/N) ~ 
0.50 for NS-BH, and (S/N) ~ 0.20 for BH-BH binaries in 
the cosmological-constant cosmology and considering the 
correlation of two Adv. LIGO detectors. If we consider 
ET in triangular configuration, the (S/N) are at least a 
factor ~ 200 greater than those obtained for Adv. LIGO. 
The signals produced in the case Chaplygin gas are always 
lower than those produced by the cosmological-constant 
cosmology. 

We have also analyzed the nature of the GW 
background produced by those compact binaries. For our 
fiducial parameters, we verify that a continuous background, 
corresponding to a duty cycle > 1, is produced by sources 
situated at cosmological distances far from z ~ 0.5 — 0.6 for 
NS-NS binaries in the cosmological-constant cosmology. For 



the Chaplygin gas cosmology, duty cycle > 1 is obtained for 
sources far from z ~ 0.6 — 0.7. 

Considering NS-BH binaries, the nature of the 
background becomes popcorn for sources far from z ~ 
0.9 — 1.0 (z ~ 1.0 — 1.5) if the dark component of the 
Universe is the cosmological constant (Chaplygin gas). On 
the other hand, BH-BH binaries are always within the shot 
noise regime. 

We also verify the characteristics of the background 
produced by a cosmological Population of stellar core- 
collapse to form black holes. In this case, signal-to- 
noise ratios within the range from ~ 0.05 to ~ 2 
could be generated depending upon the CSFR/dark-energy 
component and the efficiency of generation of GWs used. We 
verify that this Population does not behaves as a continuous 
background. However, for sources situated at cosmological 
distances z ~ 3 (z ~ 7) we verify a transition from shot 
noise to popcorn regime if the time-scale for star formation 
(r s ) is ~ 1 Gyr and if the dark-energy component is the 
cosmological constant (Chaplygin gas). 

A stochastic background of GWs, from cosmological 
origin, is expected to arise from a superposition of 
different gravitational wave sources at different redshifts. In 
particular, the formation of different objects as, for example, 
NS-NS binaries, NS-BH binaries, BH-BH systems, core- 
collapse supernovae to form black holes, among others, are 
related to the CSFR. In this way, we determine the shape of 
a pre-galactic background of GWs considering the collective 
effect of these different objects. We obtain that a stochastic 
background of GWs could be generated in the range of 
frequency 10 Hz — 1.5 kHz with a signal-to-noise ratio ~ 3 
(~ 2) for a pair of advanced LIGO interferometers and 
if the dark-energy component is the cosmological constant 
(Chaplygin gas). 

It is worth stressing that the sensitivity of the future 
third generation of detectors, as for example the Einstein 
Telescope, could be high enough to increase the expected 
values of (S/N). For example, if we consider the ET-B 
sensitivity curve in triangular configuration the gain in 
relation to advanced LIGO would be ~ 300 — 1000. Thus, 
instruments as ET could permit to explore the epoch when 
the first stars were formed in the Universe at the end of 
the so-called 'dark-ages'. In this way, the detection and 
characterization of a stochastic background of GWs, over 
a large range in frequency, could be used as a tool for the 
study of the star formation up to redshift z ~ 20. 

Recently, |Marassi, Schneider fc Ferrari ( |2009|[20Tl| also 
have analyzed stochastic backgrounds of gravitational waves 
but with a CSFR derived by |Tornatore, Fcrrara & Sch neider] 
(20071 which includes sources up to z ~ 15. Comparing the 



results of those papers with Fig. 1 of the present work, for 
the A-CDM case, we note that both CSFRs produce similar 
results up to z ~ 3. At higher redshifts (z > 3) the present 
study predicts more sources than TFS-CSFR (Tornatore, 
Ferrara & Schneider CSFR). This happens because Pereira 
& Miranda model (PM-CSFR) incorporates more baryons 
in stars than TFS-CSFR. 

The preference for using PM-CSFR comes from the 
following points: (a) The dark energy component, with 
an equation of state dependent on time, modifies both 
the expansion factor E(z) and the growth function (see 
equations |] and [§. TFS-CSFR uses the GADGET code 
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with the A-CDM model as the background cosmology. Thus, 
it would not be possible to study consistently the case 
Chaplygin gas as dark-energy fluid for the TFS-CSFR. This 
means that we would not have to compare their results with 
those obtained here for the case Chaplygin gas. 

Note that one of our results was to show that the dark- 
energy component modifies the amplitude of the CSFR; (b) 
Recently Pereira & Miranda (20111 studied four different 



CSFRs (|Pereira fc Miranda|2010||Fardal et al.|2007||Springel 
|fc Her nquist 2003; and Hopk ins &: Beacom||2006| rto derive 
the evolution of the comoving black hole mass density. Their 
results show that PM-CSFR has a good agreement with the 
quasar luminosity density up to redshift ~ 6. On the other 
hand, as PM-CSFR produces a high number of sources at 
z > 3 then there exist an important contribution of these 
objects, formed at high redshifts, to the backgrounds studied 
in the present paper. This is an intrinsic characteristic of 



the scenario used by Pereira & Miranda ( 2010 1 to derive the 
CSFR. 

As a final point, different dark-energy scenarios 
could give different signatures on the background through 
the expansion function E(z). As the gravitational wave 
background could trace the behavior of the Universe up 
to redshift ~ 20, then in principle it should be possible 
to infer if there is a temporal dependence of the dark- 
energy equation of state. That is, if w(a) ^ 0. In this 
way, not only binary systems at lower redshifts (z < 
2 — 3) working as standard siren s (e.g., |SathyaprakasE^| 



Schutz & Van Den Broeck 2010 Zhao et al. 20111 but 



also stochastic backgrounds of gravitational waves could 
contribute for a better comprehension of the physical nature 
of the dark energy. We observe that all viable dark-energy 
fluids have similar behaviour up to z ~ 1.5 where the 
main observational data are available (e.g., SNIa and baryon 
acoustic oscillations - BAO). This fact can be inferred 
from the CSFR (see Fig. 1) where for z < 2 all models 
(cosmological constant and Chaplygin gas) have similar 
evolution. Thus, a way to identify if w(a) 7^ would 
be observe the Universe at higher redshifts. In principle, 
stochastic backgrounds of gravitational waves could be such 
observable. 
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